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o : 

The recent BaBar measurements of the 77* — > ir° transition form factor show spectacular 
^ I deviation from perturbative QCD prediction for large space- like Q 2 up to 34GeV 2 . When 

, plotted against Q 2 , Q 2 F{Q 2 ) shows steady increase with Q 2 in contrast with the flat Q 2 

behavior predicted by perturbative QCD, and at 34GeV 2 is more than 50% larger than the 
QCD prediction. Stimulated by the BaBar measurements, we revisit our previous paper on 
' the cancellation of anomaly effects in high energy processes Z° — > 7r°7, e + e _ — > 7r°7 and 

JL , apply our results to the 7*7 — > ir° transition form factor measured in the e + e~ — > e + e~ir° 

process with one highly virtual photon. We find that, the transition form factor F(Q 2 ) 
behaves as (|p) X (ln(Q 2 /m 2 )) 2 and produces a striking agreement with the BaBar data for 
Q 2 F(Q 2 ) with to = 132 MeV which also reproduces very well the CLEO data at lower Q 2 . 

r- : 

(>> . PACS numbers: 11.40.Ha 12.38.Bx 13.66.Bc 

cn ; 

The 7*7 — > 7T° transition form factor at large momentum transfer Q 2 which could be measured 
m —7- 7r°7 decay, in high energy e + e — > 7r°7 or in e + e — > e + e ir° collisions [lj where one of the 
j> \ photon is highly virtual has been the subject of studies using the quark parton picture of hadrons for 

^ ■ hard exclusive processes since the earlier days of perturbative QCD. The interest in this transition 

b : 

form factor lies in the fact that it is one of the simplest quantities to compute in QCD and relatively 
easy to measure. At Q 2 = 0, it is given by the two-photon ir° decay governed by the Adler-Bell- 
Jackiw triangle chiral anomaly 2H4L wmc h gives correctly the decay rate. At large Q 2 , short- 
distance operator expansion(OPE) 3] or perturbative QCD 6-^] predicts F(Q 2 ) ~ 2f w /Q 2 (we 
use the convention /„- = 93 MeV in this paper). The earlier CLEO data 9(] give values for F{Q 2 ) up 
to Q 2 = 8GeV 2 somewhat below the perturbative QCD(pQCD) prediction, though, with a possible 
rise for Q 2 F(Q 2 ) above 2.5 GeV 2 . Recently, the BaBar Collaboration has produced measurements 
for the transition form factor from 4 to 34 GeV 2 [ljj] which show spectacular deviation from the 
perturbative QCD prediction as seen from the data for Q 2 F(Q 2 ) which rise steadily with Q 2 in 
contrast with the rather flat behavior predicted by pQCD and is more than 50% above the QCD 



prediction at 34GeV 2 . There are a" 
shown in |(], [n]] As mentioned in 10 , 



so measurements by CELLO till up to 2.5 GeV 2 which are 

"i nn 

12J , recent calculations [13, HJ] using the light-cone sum rules 



method at next-to-leading order with various forms for the pion distribution amplitude, seem to 
obtain values for the transition form factor higher than the asymptotic limit of [(], Q], but with 
very different Q 2 behavior than the BaBar data for Q 2 < 15 GeV 2 and are below the BaBar data 
for Q 2 in the range from 20 to 40 GeV 2 A more recent work |15j] seems to obtain results 

consistent with the BaBar data for Q 2 > 15 GeV 2 with a very broad pion distribution amplitude. 
Most of these calculations use the short-distance OPE of Lepage-Brodsky [{], Q] to obtain the pion 
distribution amplitude, but, in the case of the 7* — > 771"° transition form factor with one virtual 
photon with large q 2 , the short-distance expansion parameter uj = —2p- q/q 2 = 1 (p, q being the 
pion and virtual photon momentum, respectively), which is large and the calculation of Lepage 
and Brodsky for fixed but large Q 2 cannot be trusted because the corrections are important as 
mentioned in 16]. 



191 ] which 



Without further questioning the validity of the perturbative QCD prediction [171 - 
is based on the quark parton picture of the pion, one could already consider the role of chiral 
anomaly for processes involving a pion and highly virtual photons and the radiative decays of 
gauge bosons like Z° — > 7r°7 or — > 7r ± 7 decays 



201 ] . In a previous paper 



21 



22j, we showed 



that, for these processes, from the modified PCAC equation due to the Adler-Bell-Jackiw anomaly, 
the quark-parton contribution to the axial current divergence given by the triangle graph cancels 
the anomaly term resulting in the suppression of the Z° — > 7r°7 decay amplitude as well as the 
7* — > 77T° transition form factor measured in e + e~ — > 7r°7 . We also showed that for Q 2 large 
as in Z° — > 7r°7 decay or in e + e~ — > 7r°7 process with a highly virtual photon, the amplitude 
behaves like (ln(Q 2 / m 2 )) 2 / Q 2 , and a similar behavior for space-like Q 2 . The quantity Q 2 F(Q 2 ) is 
found to rise with Q 2 as (\n(Q 2 /m 2 )) 2 . This can be seen from the ln(Q 2 /m 2 )/Q 2 behavior at large 
Q 2 of the absorptive part of the triangle graph contribution to the divergence of the axial vector 
current matrix element < 7r°|3 /J yl M |7* 7 >. Hence the (ln(Q 2 /m 2 )) 2 /Q 2 behavior for the 7*7 — > ir° 
transition form factor. This is one of the main differences between the chiral anomaly approach 
and the perturbative QCD quark parton approach to the 7*7 — > tt° transition form factor. In the 
chiral anomaly approach, the triangle graph gives us the absorptive and the dispersive part, while 
the quark parton approach based on short-distance operator expansion gives us only the real part 
as given by the tree graph at the lowest order in perturbative QCD. 

Our calculation of the Adler-Bell-Jackiw triangle anomaly contribution to the 7*7 — > it tran- 
sition form factor and similar calculations later on {23], Q] are done at a time when few data on 



the transition form factors at large Q 2 are available |(], . Now that the Babar data are available 
over a large range of large momentum transfer Q 2 , it is relevant to compare data with the anomaly 
contribution, considering the fact that, in pQCD the 7*7 — > 7r° transition form factor depends on 
the pion distribution which is not known to a good accuracy at present. For this reason, in this 
paper, we apply our previous analysis of the anomaly effects in Z u — > to the 7*7 — > ir transi- 
tion form factor at space-like Q 2 . Using the generalized Goldberger-Treiman relation to relate the 
quark mass in the triangle graph to the quark-pion Yukawa coupling g, with m = g f n , as in the 



linear sigma-model 
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261 ]. we show that the BaBar data can be reproduced with m = 132 MeV 
or g = \/2, consistent with the (ln(Q 2 /m 2 )) 2 /Q 2 behavior. 

We begin by first recalling, for convenience, our derivation of the anomaly contribution to Z° — > 
7r°7 in [2l|]. Similar to the Z° — > 7r°7 decay, the 7*7 — > ir° amplitude, M =< 7r°(p)|T|7*(g)7(A;) > 
in the reaction e + e~ — > e + e~ir° has the form e pL {q)e u {k)N^ u {q,k) with: 

N^(q, k) = e 2 F{q, k)Y^, = e^q a k p . (1) 

where p is the produced pion momentum in the final state (p = q + k), F(q, k) in the following, will 
be written as F(Q 2 ), the form factor in the kinematic region of the BaBar measurement, with the 
virtual photon with momentum q space-like (q 2 = — Q 2 < 0) while the photon with momentum k is 
almost on the mass-shell (k 2 ~ 0) . As with the derivation of the two-photon 7r° decay amplitude, 
we start with the modified PCAC equation due to the Adler-Bell-Jackiw triangle anomaly in the 
presence of the electromagnetic interactions. The divergence of the axial vector current associated 
with 7T° becomes: 

8^ = Uml<P + Sr^e aM& F^F^ (2) 

with F^v the usual electromagnetic field strength tensor and S are the sum of the squares of the 
charges and colors of quark contributing to the anomaly, which takes the value S = 1/2 [2]. Taking 
the matrix element of the l.h.s of Eq. ([2]), and seperating the ir° pole term from the continuum, as 



previously shown 



21], we arrive at the expression for the 7*7 — > ir° amplitude: 



1 / „ e 2 



where R^ T (q, k) is the triangle graph (the direct coupling between the three currents) or the con- 
tinuum contribution to the axial vector current matrix element < 01^17*7 > defined as R^ UT (q, k): 

R^(q, k) = R^(q, k) - f /^^ (4) 



As shown in 0,13], gauge invariance and Bose symmetry tell us that the divergence p T R^ VT (q, k) 
is in general proportional to q 2 and k 2 and does not vanish when one or both photons are off mass- 
shell. Only when both photons are real(q 2 = 0, k 2 = 0) that p T R^ UT (q,k) is 0(p 2 ) and becomes 
negligible. One can then apply Eq. ([5]) to 7T° — >■ 77 and finds that it is given by the anomaly j3"4j. 



For highly virtual photon as in the present 7*7 — > ir° transition form factor, we wi 
[20I ]. that R^i^q^k) is given by the triangle graph. From the expression given in 2?J 



m 

assuming equal mass for u, d quarks in the triangle graph: 



1 assume, as 
we find, 



p T R flUT (q, k) = e 2 S^2mP(q,k) + ^2 (5) 



where 



(6) 



and 

D = k 2 y(l — y) + q 2 x(l — x) — 2q ■ kxy — m 2 (7) 

The quark mass m in the triangle graph is taken as a parameter to set the scale for the high energy 
limit, similar to the quark mass parameter used in the calculation of e + e~ —> qq processes at high 
energy, where the asymptotic limit is reached when Q 2 3> m 2 and the quark-parton picture is valid 
and the m 2 /Q 2 term in the cross section a(e + e~ — > hadrons) can be neglected. 
When both photons are real (q 2 = 0, k 2 = 0), from Eq. ([6]), we get: 

2mP(q,k) = -^+0(p 2 ) (8) 

which implies that the r.h.s of Eq. ([5]) is 0(p 2 ) in agreement with our previous remark that 
p T R^ UT (q,k) = 0{p 2 ). For our transition form factor with time-like virtual photon Q 2 = s, with 



s > Am 2 , we have, as given in 



2l|: 



2mP(q,k) = ^(^jK(m 2 ,s) (9) 

where 

K(m 2 , s) = ( ln |~ ~ i7r ) ' P = Vi-^/s, s>Am 2 (10) 

For space-like q, with q 2 = —Q 2 (s = —Q 2 ), with s < 0, by analytic continuation, the function 
K{m 2 , s) becomes real and is given by: 

K(m 2 ,Q 2 )= (m^j) 2 , p = y/l + Am 2 /Q 2 (11) 
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Using Eq. ([9]) for 2mP(q,k) and Eq. ([5]) for the divergence p T R^ U7 (q,k), we arrive at the final 
expression for 7*7 — > ir° amplitude: 



1 e 2 



S[^ I K(m 2 ,Q 2 )\Y^ 



(12) 



f«2v* \Q 2 ' 

We note also that the term 2mP{q, k) in Eq. ([5]) can be obtained directly from the triangle 
graph with the axial vector current replaced by the direct pion-quark vertex with the 75 Yukawa 
coupling g = m/f n as in the linear sigma model 



25 



26] and PCAC is imposed on the pion-quark 
vertex . This is an equivalent method to obtain the 7*7 — > ir° amplitude for large Q 2 without 



281 ]. The triangle graph with pion-quark 



having to go through the proof of anomaly cancellation 
vertex also gives us the term 2mP(q, k) for the 7r° — > 77 decay amplitude, as shown in Eq. 
For S = 1/2, the scalar transition form factor F(q, k) is given by: 

F(q,k) = ±-^—K(m 2 ,s) 



(13) 



f n iTT 2 S 

and for the space-like 7*7 — > ir° transition form factor, at large Q 2 m 2 , the dominant term in 
Eq. Cd is 



F(Q 2 



1 1 mV, Q 2x2 
In — 



(14) 



to 
in 



U 4vr 2 Q 2 

ae compared with the transition form factor for real photon(our normalization is the same as 



F(q 2 = 0,k 2 = 0,p 2 = 0) 



(15) 



M 2 U, 

We emphasize that our results Eqs. (|UHl5p are exact calculations and the result for the transition 
form factor given in Eq. (|14p is valid for large Q 2 ( Q 2 3> m 2 ), including its Q 2 — > 00 limit. 

This result, the (ln[Q 2 /m 2 )) 2 rise for Q 2 F(Q 2 ) which has been obtained in [2l|], are obtained 



later in 



Rcf. 



23] obtains the anomaly contribution from the PCAC equation for the matrix 



element of the axial vector current between the photons, while Ref. [24J equates the divergence of 
the axial vector current two-photon matrix element with its pion pole terms. In our calculation, 
the quark mass parameter is taken as a dynamical (constituent) quark mass. The reason why we 



lake const il uenl quark mass is also emphasized in |2^,Q. In fact as argued in Rcf. 23 . because 
of the presence of the Nambu-Goldstone pion pole term in the PCAC equation, the divergence 
of the axial vector current should also be non-vanishing in the limit of vanishing current quark 
mass, to cancel the longitudinal term generated by the pion pole term, as in the derivation of the 
Goldberger-Treiman relation for the pion-nucleon coupling constant 
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As shown in Fig. (pQ), our prediction for the quantity Q 2 F(Q 2 ) for m = 132 MeV fits very well the 
CLEO and BaBar data. The agreement with the BaBar data is striking, as our predicted transition 
form factor depends on only one parameter, the effective mass for quark in the triangle graph . We 



13, 



151 ] seem 



note that recent works with various models for the pion distribution mentioned above 
unable to obtain the rise of Q 2 F(Q 2 ) for Q 2 > 20 GeV 2 as the BaBar measurements. Our value for 
the effective quark mass is consistent with the high energy behavior of the e + e~ — > hadrons cross 
section for which the ratio R = a{e + e~ —> hadrons) /a(e + e~ —> is constant for s above a 

few GeV 2 . 
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FIG. 1: Chiral anomaly prediction(solid line) for Q 2 F(Q 2 ) compared with the BaBar and CLEO measured 
values and the large Q 2 pQCD prediction (horizontal dash line) of [6| 



In conclusion, we have shown that chiral anomaly effects produce a (m 2 / Q 2 ){\n(Q 2 / m 2 )) 2 be- 
havior for the 77* — > 7r° transition form factor at Q 2 S> m 2 in contrast with the 2f w /Q 2 behavior 
given by perturbative QCD. It is remarkable that our simple expression for the transition form 
factor is able to explain the CLEO and BaBar data. Our prediction for the rise of Q 2 F(Q 2 ) at 
higher space-like Q 2 could be confirmed with further measurements. Similar behavior is expected 
for the time-like transition form factor and could be seen in e + e~ —> vr°7, as in e + e~ — > 77^7 for 
which the time-like transition form factor has been measured by CLEO at Q 2 = 14.2 GeV 2 [30j and 
is found to be close to the corresponding space-like values 12|, and by BaBar at Q 2 = 112 GeV 2 
3l| . So, if chiral anomaly is indeed the cause of the rise like (ln(Q 2 /m 2 )) 2 of Q 2 F(Q 2 ), then, the 
pion might very well, for 77* — > 7r° process at high energies, behave like a Nambu-Goldstone boson, 



7 



like the Adler zero in low-energy tttt scatterings and in ip' — > J/ i/jittt decay which are obtained from 
chiral symmetry constrains. 

Note added. After the completion of this paper, we were informed of the papers by Dorokhov 
331 ] in which similar results are obtained using the triangle graph with pion-quark coupling 
given in j^ . 
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